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Abstract 

A single - phase melting moving boundary problem is formulated after a 
change of variable, as an parabolic variational inequality. 

Introduction 

The term moving boundary problems (MBP'S) is commonly used when the 
boundary is associated with time dependent problems and the boundary of the domain 
is not known in advanced but has to be determined as a function of time and space. 

Heat conduction or diffusion type problems with phase change constitute a 
large class of (MBP'S) usually called Stefan problems. Moving boundary problems 
have received much attention in recent year due to their practical importance in 
engineering and science [1], [2], [3]. 

The purpose of this paper is to obtain a variational inequality formulation of a 
single-phase melting problem. 

Description of the problem 

The problem concerns heat transfer in an ice-water medium occupying the 
region 0 < X < 1 . At any time t, the water that under going phase change, is 
contained in the region and the rest of the region outside it, is occupied 

by ice, as in [4]. The temperature of the water is assumed to be equal to zero which is 
also the critical temperature of phase change. 

Mathematical Formulation of the problem 

The problem is expressed in non-dimensional form. If u(x,t) denotes the 
temperature distribution at distance x and time t, s(t) is the position of the interface at 
time t, then the problem is to find the pair of unknowns u(x,t) and s(t). The problem 
may be defined mathematically by the following equation and conditions [5]. 
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the boundary conditions are given by :- 

x=0 , t > 0 (2) 

x>0 , t=0 (3) 



u=0 
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S(0) = 1 



.(5) 



Where A is "Latent heat" and (4) is known as the "Stefan condition". 



Variational Inequality Formulation of the problem 

Another possibility is to formulate moving boundary problems in terms of 
certain inequalities, either of a differential or variational nature. 

In either case the expressions refer to a fixed domain and explicit use of the 
Stefan or other interface condition is avoided. A variational formulation is only 
possible after some suitable transformation has removed the discontinuity. Such a 
transformation was first proposed by Duvaut (1973) for moving boundaries prompted 
by an idea of Baiocchi (1971) for stationary, free boundaries. Fremoud (1974, 1980) 
called the new variable , the freezing index. 

For the single-phase melting problem defined by equations (1-5), the Duvaut 
transformation, or the Fremond freezing index, is 

w[x,t)- ju{x,r)dT ,0<x<s(t) 

t(x) 

(6) 

w(x,t)=0 ,s(t)<x<l 

where was have introduced the fixed, finite domain 0 < X < 1 , and where 
t — £{x) define the time when the ice at point x changes to water, i.e. £{x) is the 

time when the melting boundary is at x. It is the inverse of s(t) so that £ 1 (?) = s(t^). 

dw 

To show that w(x,t) and its first derivative are continuous on the melting 

dx 

boundary and throughout the range 0 < X < 1 , we have form (6) that:- 

— = f — (x,r)dT- — u(x,£(x)) 0<x<s(t) 

dx ( ( x ) dx dx 

— = 0 s(t)<x<l (7) 

d x 

since the temperature u is zero on the interface x=s(t) , i.e, u(x,^(x))=0, then (7), 
becomes :- 

— = J — (x,r)dr ,0<x<s(t) 

dx ( ( X ) dx 

(8) 



dw 

— = 0 ,s(t)<x<l 

dx 

dw 

and form (8) we know that = U on x=s(t). 

dx 

Furthermore , 

d 2 w f d 2 u , didu . . . ~ / \ (n\ 

— = j —dT- — — (x,£(x)), 0<x<s(t) (9) 

ox f{x) ox Ox Ox 
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Noting that — = 

dx di 
dt 

So that the Stefan condition (4) implies 

du dl 

= A , then the equation (9) with use of (1) becomes 

dx dx 
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Thus the transformed equation is :- 

$ w d 2 w n ^ 

= - + A ,0<X<s(t) 

dt dx 2 



(10) 



w=0 , s(t) <x<l 

with w(x,t) > 0 in 0 < X < s(t) . 

Thus w(x,t) and its first derivative are continuous throughout the domain 0 < X < 1 , 
then we have the inequalities :- 

dw d 2 w 



dt dx 2 
dw d 2 w 



and 
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-A>0 

-A\w = 0 



,w>0,0<x<L 



,0<x<l- 
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Avariational form of (11) and (12) is derived by multiplying the left-hand side of the 
first of (110 by (v-w) using the test functions v which are non-negative in 0<x<l with 

dv 

v=0 on x=l, and for which v and — are square integrable in 0<x<l, i.e 

dx 

jv 2 dx and J [ — | dx are bounded 



o 

Thus 

dw d 2 w 



dt dx 2 
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-X l(v - w)dx 



(v - w) 
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+ J — — (v - w)dx -\a{v - w)dx 

o o dx dx o 
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Remembering that v=w=0, x=l, and = 0, x=0, together with use of (12), we 

dx 

obtain after slight rearrangement :- 

dw d 2 w 



Aj(v- w)dx + J 



dt dx 2 



-I 



vdx 



But the final term on the right-hand side of this equation is non-negative since v > 0, 

and the first of (1 1) is to be satisfied. 

We therefore have the variational inequality 

J — (v - w)dx + j — — (v - w)dx > Aj(v - w)dx (13) 

o dt o ux o 

Define the inner products (.,.) and a(.,.) with reference to a domain D in two space 
dimensions as follows :- 

a(u, v - u) = j" j" VwV(v - u)dx dy {l4) 

D 

= ll{ U A V X - U X ) + U y( V y- U y )} dx 
D 

and (f,v-u) = JJ f(v -u)dx dy (l5) 

D 

using the nomenclature of (14) and (15) , the inequality formulation of the melting 
problem (13), equivalent to the differential forms (11) and (12), becomes 

(w t ,v- w)+ a(w,v - w)> (A,v- w) (16) 
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